Experimental Validation of the Largest Calculated Isospin-Symmetry-Breaking Effect 

in a Superallowed Fermi Decay 
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A precision measurement of the 7 yields following the /3 decay of 32 CI has determined its isobaric 
analogue branch to be (22.471q' 2 ,|,)%. Since it is an almost pure Fermi decay, we can also determine 
the amount of isospin-symmetry breaking in this superallowed transition. We find a very large 
value, 5c = 5.3(9)%, in agreement with a shell-model calculation. This result sets a benchmark 
for isospin-symmetry-breaking calculations and lends support for similarly-calculated, yet smaller, 
corrections that are currently applied to + — > + transitions for tests of the Standard Model. 



PACS numbers: 24.80.+y, 23.40.Bw, 29.30.Kv 

Precisely measured ft values of J 71 " = + — > + j3 de- 
cays of isospin T = 1 nuclei are used to set stringent 
limits on scalar and right-handed interactions, verify 
the conserved vector current (CVC) hypothesis at the 
~ 10 -4 level, and provide the most precise measurement 
of V u d, the up-down clement of the Cabibbo-Kobayashi- 
Maskawa quark- mixing matrix 0, Q. For these pur- 
poses, the vector coupling constant is extracted from the 
experimental ft values of these nuclei, after correcting 
for isospin-symmetry-breaking and radiative effects. To 
date, the ft values of 13 such nuclei have been experimen- 
tally determined to a precision of < 0.3% which sets one 
of the most demanding tests of the Standard Model 0- 
These experimental advances have placed the theoreti- 
cally calculated corrections under intense scrutiny in re- 
cent years. Emphasis has been placed on the nuclear- 
structure-dependent isospin symmetry breaking (ISB) 
corrections, denoted by 5 c an d defined by the equa- 

tion |A/ F | 2 = \M \ 2 (l-S c )- Here M F is the Fermi matrix 
element for the transition and M is its value in the limit 
of strict isospin symmetry, which is broken by Coulomb 
and charge-dependent nuclear forces. The 13 cases just 
mentioned are all T = 1 — > 1 transitions in A = An + 2 
nuclei with corrections that are small and of order 1%. 

A discriminating test of these calculations is realized by 
investigating cases where the correction is much larger. 
Up to now, though, there have been no nuclei studied 
where 5c is larger than ~ 2% @, 0, @. In this Letter, 
we focus on the /3 decay of 1 + ,T = 1 32 C1 as a test 
of isospin-mixing calculations. Its Fermi decay branch 
feeds the analogue 1 + , T = 1 state in 32 S, whose position 
in the spectrum at 7002-kcV excitation is very close to 
a known 1+,T = state at 7190 keV Q. This greatly 
enhances the size of the isospin-breaking correction. Our 
calculation of the ISB effect for this A = An nucleus 
is 5c = 4.6(5)%, a value significantly larger than those 
found in any of the A = An + 2 nuclei. Thus the case 



of 32 CI provides a unique opportunity to test isospin- 
symmetry-breaking calculations where the correction is 
relatively very large. 

The experiment was performed at the Cyclotron In- 
stitute at Texas A&M University. Details of this exper- 
iment will appear in a separate paper Q. Briefly, we 
produced 32 CI via the inverse- kinematic transfer reaction 
1 H( 32 S, n) 32 Cl using a LN 2 -cooled, H 2 gas target with a 
400 nA 32 S primary beam at 24.8 MeV/nucleon. The 
reaction products were spatially separated by the Mo- 
mentum Achromatic Recoil Separator [||, resulting in a 
91% pure, 20 MeV/nucleon 32 C1 beam with an intensity 
of ~ 2 x 10 5 ions/s. The beam was implanted and col- 
lected in an aluminized-Mylar tape for 0.8 s before a fast 
tape-transport system moved the activity to a shielded 
counting station 90 cm away. Data for f3 — 7 coincident 
events were acquired using a 1.5 inch diameter, 1 mm- 
thick scintillator and a 70% HPGc detector. Count times 
were for 1, 2 and 4 sec (76%, 13% and 11% of the data 
respectively). The scintillator was placed 0.5 cm from 
the activity, detecting > 40 keV positrons with « 32% 
efficiency. On the opposite side of the tape the HPGc 
was placed a large distance away (15.1 cm) to reduce the 
effects of coincidence summing of the 7 rays. The cycle of 
collecting, transporting and measuring the 32 CI activity 
was repeated continuously throughout the experiment. 

Critical to the success of this experiment was the ex- 
tremely precise photopcak efficiency calibration of our 
HPGe detector. As described in detail in Ref. [l(|, using 
a combination of measurements and Monte Carlo calcu- 
lations with the code Cyltran 11 



tcrmincd to ±0.2% from £L = 50 - 



, the efficiency is dc- 
1400 keV and ±0.4% 



from 1.4 — 3.5 MeV. The energy range of the HPGe in the 
present work, however, extends up to 7.35 MeV; so we ex- 
tended the efficiency curve above 3.5 MeV using the same 
Cyltran code used in Ref. [13]. To estimate uncertain- 
ties in this extrapolation, we performed an independent 
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FIG. 1. Top: Percent differences in the absolute effi- 
ciency from Ref. [Tol ] (solid line) and our adopted uncer- 
tainties (shaded region). This adopted efficiency curve and 
its CYLTRAN-calculated extrapolation (dashed line) are com- 
pared to a Penelope simulation (dotted line). Bottom: 7- 
energy spectrum of the 11.2 x 10 6 /3-coincident events ob- 
served in this work. Aside from the isolated 30 S contaminant 
peak at 677 keV, labeled peaks are associated with the de- 
cay of 32 C1 (* and ** indicate their single- and double-escape 
peaks). Boxed values indicate transitions from the 7002-keV 
isobaric analogue state. 



calculation using the Monte Carlo code Penelope [12| . 
The difference between the two efficiency curves is shown 
in the top panel of Fig.Q] Over the range of measured val- 
ues, Penelope reproduces well the experimentally deter- 
mined efficiency, mostly within the ±(0.2 — 0.4)% uncer- 
tainty range from Ref. [10(; above 3.5 McV, the difference 
between the two extrapolations is contained within our 
adopted uncertainties ranges of ±1% from 3.5 — 5 MeV 
and ±5% above 5 MeV. The bottom panel of Fig. Q] 
shows a plot of the observed 7 spectrum in coincidence 
with a j3 signal in the scintillator, where nearly every 
observed peak is associated with the decay of 32 CI. 

Once the peak areas were obtained, we used the pre- 
cisely known efficiencies to convert them into relative 
yields of 7-rays. We then fit the (3 and 7 branches to re- 
produce these yields. Our measurement has found 3 new 
P branches, 22 new 7 lines, placed limits on 10 potential 
7 transitions, and improved the precision of the branches 
and yields reported previously [13j by about an order of 
magnitude [8j. The 12/3 transitions we observe and the 
known ground state branch of (l.Ojlo 2)% from Armini et 
al. [IH represent almost all of the (3 yield; however, there 
is potentially still a large number of weak (3 transitions 
which, though too weak to be seen individually, may sum 
up to a total /? strength that is non- negligible. This "Pan- 
demonium effect" [IH was recently raised [HI in the con- 
text of superallowed f3 decay in p, /-shell nuclei. Here we 
follow the approach advocated there: to compute these 
very weak, unobserved j3 branches using a shell-model 



calculation and include this predicted strength as a small 
correction in the analysis. We take the model space to 
be the full s,d shell and use the charge- independent ef- 
fective interaction of Wildcnthal USD 17J , as well as the 
more recent USDA and USDB updates 18| of Brown and 
Richter pi]. 

Our analysis of the branches and yields includes a to- 
tal of 51 excited states in 32 S. In addition to the 12 
(3 branches observed, our shell model calculation iden- 
tifies approximately 30 additional weak (3 transitions to 
states whose excitation energy lies between 7.485 and 
« 11.8 McV. Though none of these individually has a 
/3-transition strength greater than 0.3%, the summed (3 
strength of all of them is 0.60(10)%, where the uncer- 
tainty is a result of the different interactions used in the 
shell model. We include these weak (3 strengths and the 
de-excitation 7 rays predicted by the shell model in our 
overall analysis to account for the small Pandemonium 
effect. We do not separately include the a-particle and 
proton-emitting states reported by Honkanen et al. [l9| 
because their summed f3 strength of 0.080(5)% is signif- 
icantly smaller than-and is no doubt already included 
in-the missing strength predicted by the shell model. 

We find the (3 branch to the isobaric analogue state 
(IAS) at 7002 keV is R = (22.47 ±0.13^;^)%. The first 
uncertainty is statistical and the second is dominated by 
two sources of systematic uncertainty: fo'ol% fr° m the 
(1-0^q'5)% ground state branch reported by Armini et at, 
and ±0.10% from the photopcak efficiency of the HPGe 
detector. 

To derive the experimental ft value, we obtain the 
partial half-life, t, from 
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Pec) = 1.327(13) s, 



(1) 



where the 32 C1 half-life is t 1/2 = 298(1) ms [14J, R is 
the superallowed branching ratio quoted above, and the 
small electron-capture fraction is calculated to be Pec = 
0.071%. We use the shell model to compute the shape 
correction function C(W) (as described in the appendix 
of Ref. [l[ ) when defining the statistical rate function 
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P W(W ~W) 2 F(Z,W)C{W)dW, (2) 



where W = E e /m e is the total energy of the positron 
in electron rest-mass units, p = (W 2 — l) 1 / 2 is its mo- 
mentum, Z is the charge of the daughter nucleus, and 
F{Z 1 W) is the Fermi function. The end-point energy, 
Wq, is determined usin g — 26015.535(2) keV for the mass 



32 ci. 



excess of 32 S from Ref. [20j, and we average Refs. 
to get —13334.60(57) keV for the mass excess of 
Combined, the decay energy is Qec = 12680.9(6) keV. 
This gives / = 2411.6 ± 2.3 ± 0.3 for the phase-space fac- 
tor, where the first uncertainty is from the Qec value and 
the second is from the shell- model calculation of C(W). 
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Thus the experimental ft value for decay to the IAS is 
ft = 3200(30) s, where the precision is dominated by the 
±0.9% uncertainty in the branch to the IAS. 

We now deduce an experimental value for the isospin- 
mixing parameter, 8c, of the decay to the isobaric- 
analogue state from the measured ft value. Before this 
can be done, however, we need some way of separating 
out the Gamow- Teller component of this mixed 1 + — > l + 
transition so that we can analyze the Fermi component 
alone. Fortunately, the USD, USDA and USDB shell-model 
calculations (described later) all predict that the Gamow- 
Tellcr matrix element is « 0.1% of the Fermi matrix ele- 
ment, and so is negligibly small for this transition. Thus 
we proceed to analyze this transition as if it were a pure- 
Fermi type via the equation 



ft(i + s' R )(i 

Here K/{hcf = 



JNS 



Sc) = 



B(F) + B(GT) 



(3) 



27r 3 fi,ln2/(m e c 2 ) 5 is a constant and 



Gv is the vector coupling constant characterizing the 
strength of the vector weak interaction. The numera- 
tor in Eq. Q may be evaluated using the precision work 
on + — > 0+ superallowed transitions: K/Gy(l + A^) = 
2 {Ft _>0 ), where the average of the 13 most precisely- 
measured cases yields (J"i 0+ ^ 0+ ) = 3071.81(83) s @- The 
radiative correction has been split into three pieces: (a) a 
nucleus-independent term, A^ , is included in { J-t° ) ; 
(b) a trivially nucleus-dependent term, S' R , is calculated 
to be 1.421(32)%; and (c) a second nuclear-structure- 
dependent term, <Sns- is determined to be —0.15(2)% 
in a shell-model calculation following the procedures in 
Ref. [H]. Finally B(F) and B(GT) are the squares of the 
Fermi and Gamow- Teller matrix elements. 

In the isospin-symmetry limit, B(F) = |Mo| 2 = 2 for 
T = 1 transitions. For B(GT), we take the three theo- 
retical values from the shell model using the USD, USDA 
and USDB effective interactions, average them and assign 
an uncertainty which spans the three calculated values: 
B(GT) = (l.Stl 3 7 ) x 10 ~ 3 - This is negligibly small com- 
pared to the dominant Fermi strength. On rearranging 
Eq. ([3]), we obtain 



<^ xp = 1 



Sns - 
5.3(9)%. 



2(Tt ' 



ft(l + S' R )[B(F)+B(GT)] 



(4) 



This represents a very substantial isospin-symmetry- 
breaking term, the largest ever determined in a super- 
allowed Fermi transition. As Fig. [5] shows, it is an order 
of magnitude larger than the typical correction applied 
to the + — >■ + pure Fermi decays, and nearly 3x larger 
than the biggest of these cases, 74 Rb. Thus this provides 
a strong benchmark with which to compare the variety 
of theoretical methods and models proposed to calculate 
ISB in nuclei 
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FIG. 2. Our determination of the isospin-breaking correction 
for 32 CI (filled circle), and calculations for 32 CI as well as other 
superallowed transitions (open points), with the Sci and 5c2 
components shown separately. The measurement and predic- 
tion for 32 C1, particularly the Sci component, is significantly 
larger than in any of the + — > + transitions. 



In what follows we present a shell-model calculation of 
dc to compare to the above result following the proce- 
dures developed by Towner and Hardy Q ■ The technique 
is to introduce Coulomb and other charge-dependent 
terms into the shell-model Hamiltonian. However, be- 
cause the Coulomb force is long range, the shell-model 
space has to be very large indeed to include all the po- 
tential states that the Coulomb interaction might con- 
nect. Currently this is not a practical proposition. To 
proceed, Towner and Hardy divide 5c into two parts: 
5c = 5ci +5c2, where Sci arises from configuration mix- 
ing between states of the same spin in a shell-model cal- 
culation using a restricted basis (in this case the full s, d 
shell), while Sc2 separately encompasses mixing beyond 
this model space. 

Starting with Sci, we perform a shell- model calcu- 
lation in the truncated Ofouj model space of the s, d- 
shell orbitals. Charge-dependent terms are added to 
the charge-independent Hamiltonians of USD, USDA and 
USDB. The strengths of these charge-dependent terms 
are adjusted to reproduce the b = —5.4872(35) MeV 
and c = 0.1953(37) MeV ^ coefficients of the isobaric 
multiplct mass equation as applied to the 1 + ,T = 1 
triplet of states in A = 32, the states involved in the (3- 
transition under study. As already mentioned, the bulk 
of the isospin mixing in the IAS occurs with the neigh- 
boring 1 + ,T = state. In the limit of two-state mix- 
ing, perturbation theory implies that Sci oc 1/(A_E) 2 , 
where AE is the energy separation of the analog and 
non-analog 1 + states. Experimentally, it is known to 
be 188.2 ± 1.2 keV @, H (compared to the much larger 
2-4 MeV of most 0+-> 0+ transitions 9). The shell 
model calculates this separation to be 184 keV with USD, 
248 keV with USDA and 387 keV with USDB interactions. 
We avoid the large uncertainties this would impose on our 
calculation by following the Towner-Hardy recommenda- 
tion Q of scaling the calculated Sci value by a factor 
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of (Ai?) 2 hco /(A_E) 2 xp , the ratio of the square of the en- 
ergy separation of the 1 + states in the model calculation 
to that known experimentally. Following this procedure, 
the 8ci values obtained in the three shell-model calcu- 
lations are reasonably consistent: Sci = 3.73% for USD, 
3.32% for USDA, and 4.19% for USDB. We average these 
three results and assign an uncertainty equal to half the 
spread between them to arrive at dci = 3.75(45)%. As 
Fig. [2] shows, this is over an order of magnitude larger 
than Sci calculated for any of the thirteen + — > + tran- 
sitions used to determine V u d- 

For the calculation of 5c2 we consider mixing with 
states outside the Ofiw shell-model space. The princi- 
pal mixing is with states that have one more radial node. 
Such mixing effectively changes the radial function of the 
proton involved in the (3 decay relative to that of the 
neutron. The practical calculation, therefore, involves 
computing radial overlap integrals with modeled proton 
and neutron radial functions. Details of how this is done 
are given in Rcf. Q. The radial functions are taken 
to be eigenfunctions of a Saxon- Woods potential whose 
strength is adjusted so that the asymptotic form of the 
radial function has the correct dependence on the sep- 
aration energy. The initial and final A-body states are 
expanded in a complete set of (A — l)-parent states. The 
separation energies are the energy differences between 
the A-body state and the (A — l)-body parent states. A 
shell-model calculation is required to give the spectrum 
of parent states and the spectroscopic amplitudes of the 
expansion. For the three USD interactions, we compute 
S C 2 = 0.827% for USD and 0.865% for both USDA and 
USDB. Our adopted value is S C 2 = 0.85(3)%. The un- 
certainty, calculated in the same manner as described in 
Ref. [H, represents the range of results for the USD in- 
teractions, the different methodologies considered in ad- 
justing the strength of the Saxon- Woods potential, and 
the uncertainty in the Saxon- Woods radius parameter as 
fitted to the experimental charge radius of 32 S. 

Combining our adopted shell-model calculations, 
5 C i = 3.75(45)% and 5 C2 = 0.85(3)%, we find <^ hcor = 
4.6(5)%, which agrees with the experimentally deter- 
mined 5.3(9)% of Eq. (|4|) within stated uncertainties. 
The agreement between theory and experiment in this 
case where 8c is so large represents a very impor- 
tant validation of the theoretical procedures outlined 
here to calculate the ISB effects in nuclei. In par- 
ticular, for (shell-model) calculations which separate 
configuration-mixing and radial-overlap components, this 
<5ci-dominated result provides an especially sensitive 
benchmark for the approximations used when calculat- 
ing configuration-mixing contributions to the total ISB 
effect in superallowed + — > + decays. 

In conclusion, we have measured relative 7-ray intensi- 
ties and characterized the superallowed (3 branch for the 
decay of 32 CI. The isospin-symmetry-breaking correc- 
tion, 5c , of this almost pure Fermi transition is consider- 



ably larger than the typical values found in other super- 
allowed decays and can therefore be used as a stringent 
test of theoretical procedures to calculate these isospin- 
brcaking effects. Agreement with a shell-model calcula- 
tion validates the approach currently applied Q to the 13 
most precisely measured ft values of + — > + transitions 
and adds confidence to the small corrections applied to 
them for Standard Model tests. 
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